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Itl  AdlTW  C  T 


The  detection  of  spread-spectrum  multiple-access  signals  relayed  through  a  hard- 
limiting  repeater  is  of  considerable  oractical  interest  in  satellite  communication. 
Analysis  of  the  nroblem  is  difficult  because  of  the  nonlinear  characteristic  of  the 
limiter.  The  approach  generally  taken  is  to  calculate  the  receiver  output  signab^-tb- 
noise  ratio  (SNR)  and  then  to  use  results  relating  error  probability  and  SNR,  which 
are  strictly  valid  only  for  a  linear  channel. 

Thie^study  provides  a  method  for  analyzing  accurately  the  detection  of  a  hard- 
limited  PSK  signal  buried  in  noise  Laving  non-Gaussian  components,  by  using  a  cor¬ 
relation  detector.  The  problem  is  of  both  theoretical  and  practical  interest.  The 
probability  of  error  is  shown  to  be  directly  calculable  from  the  characteristic 
function  of  the  receiver  output,  without  the  need  of  obtaining  the  probability  den¬ 
sity  function,  which  generally  requires  a  difficult  inverse  transformation  of  the 
characteristic  function.  This  approach  might  be  of  some  theoretical  interest  in 
solving  related  problems.  v  ( 

The  detection  of  a  biphaAe-modulated  spread-spectrum  signal  relayed  through  a 
hard-limiting  repeater  is  analyzed  without  resorting  to  heuristic  SNR  arguments.  An 
analytical  expression  for  the  error  probability  at  the  receiver  output  is  obtained 
by  considering  the  noise  present  at  both  the  satellite  repeater  input  ard  the  ground 
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13.  Abstract  (continued) 

receiver  input.  Numerical  evaluation  of  the  erroJ 
is  presented  for  various  values  of  processing  gair 
link  input  noise  levels. 

The  detection  of  a  binary  PSK  signal  after  ti 
limiter  is  also  analyzed.  The  error  rate  in  this 
form  of  a  series  containing  either  confluent  hypei 
modified  Bessel  functions.  Numerical  results  for 
case  and  the  PSK  case  are  compared  with  those  for 
assess  the  system  performance  degradation  caused  t 
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The  detection  ot  spread-spectrum  multiple-access  signals  rcla>ed 
through  a  hard-limiting  repeater  is  ol  considerable  practical  interest 
in  satellite  communication.  Analysis  of  the  problem  is  dil'ficuit  be¬ 
cause  of  the  nonlinear  characteristic  of  the  limiter.  The  approach 
generally  taken  is  to  calculate  the  receiver  output  signal- to-noise 
’’atio  (SNR)  and  then  to  use  results  relating  error  probability  and  SNR, 
which  are  strictly  valid  only  for  a  linear  channel. 

This  study  provides  a  method  for  analyzing  accurately  the  detection 
of  a  hard-limited  PSK  signal  buried  in  noise  having  non-Gaussian  com¬ 
ponents,  by  using  a  correlation  detector.  The  problem  is  of  both  theoret¬ 
ical  and  practical  interest.  The  probability  of  error  is  shown  to  be 
directly  calculable  from  the  characteristic  function  of  the  receiver 
output,  without  the  need  of  obtaining  the  probability  density  function, 
which  generally  requires  a  difficult  inverse  transformation  of  the 
characteristic  function.  This  approach  might  be  of  some  theoretical 
interest  in  solving  related  problems. 

The  detection  of  a  biphase-modulated  spread-spectrum  signal  relayed 
through  a  hard-limiting  repeater  is  analyzed  without  resorting  to  heuristic 
SNR  arguments.  An  analytical  expression  for  the  error  probability  at  the 
receiver  output  is  obtained  by  considering  the  noise  present  at  both  the 
satellite  repeater  input  and  t.he  ground  receiver  input.  Numerical  evalua¬ 
tion  of  the  error  probability  expression  is  presented  for  various  values 
of  processing  gain  and  up-link  and  down-link  input  noise  levels. 

The  detection  of  a  binary  PSK  signal  after  transmission  through  a 
hard  limiter  is  also  analyzed.  The  error  rate  in  this  case  is  obtained 
in  the  form  of  a  series  containing  either  confluent  hypergeometric  func¬ 
tions  or  modified  Bessel  functions.  Numerical  results  for  both  the  spread- 
spectrum  case  and  the  PSK  case  are  compared  with  those  for  a  linear  re¬ 
peater,  to  assess  the  system  performance  degradation  caused  by  limiting. 
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EVALUATION 

The  significance  of  this  report  is  that  it  derives  the  expression 
for  the  detected  error  probability  for  a  biphase  modulated  signal  trans¬ 
mitted  through  a  bandpass  hard  limiting  channel.  Tnus,  RADC  engineers 
can  obtain  the  exact  performance  of  such  a  communications  relay  channel, 
or  any  receiver  containing  &  hard  limiter.  Clearly  evident  from  this 
report  is  the  parametric  range  over  which  the  usual  linearizing  assumption 
is  valid.  The  analysis  set  forth  therein  provides  a  framework  for  exact 
future  analysis  of  mu-linear  channels  such  as  those  used  with  mo.  e 
complicated  signal  structures,  multiple  access,  and  soft  limiting 
or  other  non-linear  channels. 

«jbHH  jy PATTI,  Project  Engineer 
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I  INTRODUCTION 

The  performance  of  c  gital  communication  systems  is  usually  charac¬ 
terized  by  the  probabilil  of  error  in  the  detection  of  the  transmitted 
signal  at  the  receiver.  hen  the  system  contains  a  nonlinearity,  such 
as  a  limiter,  it  is  diffj  ult  to  evaluate  the  effects  of  the  nonlinear 
device  on  the  system  perl  rmance.  The  approach  generally  taken  is  to 
calculate  the  receiver  oi  put  signal-to-noise  ratio  (SNR),  and  then  to 
use  results  relating  the  rror  probability  and  the  SNR,  which  are  strictly 
valid  only  for  a  linear  <  annel.  A  comprehensive  investigation  of  the 
effect  of  hard  limiting  c  signal  detectability  for  a  system  consisting 
of  a  limiter,  a  narrowbar  filter,  and  an  envelope  detector  in  cascade 


is  given  in  Refs.  1  and  2.  Aein  considered  the  channel  structure  in 
which  the  narrowband  filter  and  the  envelope  detector  have  been  replaced 
by  a  correlation  detector,  and  he  analyzed  the  effect  of  limiting  on  the 
probability  of  error  in  the  detection  of  a  constant-envelope,  phase- 
coded  spread-spectrum  signal.  Aein  derived  the  expression  for  the  bit- 
error  probability  under  the  following  assumptions: 

(1)  The  SNR  at  the  limiter  input  is  small. 

(2)  The  processing  gain  (TW  product)  of  the  system  is  large. 

He  showed  that  under  the  above  limitations  the  amplitude  distribution  of 
the  interference  (the  sum  of  down-link  and  up-link  retransmitted  noise) 
at  the  receiver  output  is,  to  a  good  approximation,  Gaussian  and  that 
the  limiter  can  be  regarded  as  a  quasi-linear  device  that  degrades  the 
system  performance  by  tt/4  or  -1.05  dB. 


References  are  listed  at  the  end  of  the  report. 
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Tins  study  derives  an  expression  for  the  error  rate  without  re¬ 
sorting  to  the  simplifying  assumptions  made  by  Aein.  Consequently,  the 
bit-error  probability  expression  derived  here  is  valid  for  a  considerably 
larger  range  of  values  of  SNR  and  processing  gain.  The  analysis  shows 
that  the  amplitude  distribution  of  the  interference  also  has  non-Gaussian 
components,  which  must  be  considered  when  dealing  with  low  error  rates. 

The  expression  for  the  error  probability  is  given  by  a  series  whose 
leading  term  is  identical  with  Aein's  result  when  the  SNR  at  the  limiter 
input  is  small  and  the  processing  gain  is  large. 

The  specific  model  fo-'  the  communication  system  considered  in  this 
study  is  shown  in  Figure  1.  The  input  to  the  limiter  consists  of  a  single 
biphase-modulated,  constant-envelope,  phase-coded,  spread-spectrum  signal 
and  a  band  of  zero-mean,  stationary  Gaussian  noiso.  The  bandpass  filter 
preceding  the  limiter  is  assumed  to  be  wide  enough  to  pass  the  signal 
with  negligible  distortion  and  to  limit  the  input  noise  to  a  narrow 
bandwidth  that  is  small  compared  to  the  center  frequency  of  the  filter. 

The  limiter  has  a  hard-limiting  characteristic  that  limits  its  output  to 
either  1  or  -1.  It  is  followed  by  an  ideal  zonal  bandpass  filter  that 
confines  the  limiter  output  spectrum  essentially  only  to  the  fundamental 
band  of  the  signal.  After  passing  through  the  satellite  repeater,  the 
signal  is  transmitted  to  the  receiver,  and  independent  thermal  noise  is 
added  to  it  on  the  down  link.  The  receiver  processes  the  composite  sig¬ 
nal,  extracting  the  information-bearing  signal  through  a  correlation 
operation  with  a  locally  generated  replica  of  the  transmitted  code  of 
the  desired  signal.  It  is  assumed  that  bo'-h  chip  and  bit  synchronization 
at  the  receiver  has  been  achieved  and  is  maintained. 

The  applicability  of  the  system  model  of  Figure  1  is  much  wider  than 
the  objective  of  this  study,  which  is  detection  of  a  single  biphase-coded 
spread-spectrum  signal  at  the  receiver  after  transmission  through  a  hard- 
limited  satellite  repeater.  The  model  is  equally  valid  for  the  p.ialysis 
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of  the  error  rate  of  a  hard-ltraited  PSK  signal.  In  this  case,  the  re¬ 
ceiver  reference  signal  is  synchronized  in  both  phase  and  frequency  with 
the  desired  transmitted  signal,  for  successful  demodulation  of  the  signal 
at  the  receiver.  An  analytical  expression  for  the  error  rate  in  the  de¬ 
tection  of  a  single  PSK  signal  has  been  obtained  in  the  form  of  a  series 
containing  either  confluent  hypergeometric  functions  or  modified  Bessel 
functions.  The  expression  is  valid  for  arbitrary  values  of  up-link  and 
down- link  SNRs. 


In  the  absence  of  down-link  noise,  the  model  represents  a  system 
with  a  limiting  front  end  in  a  receiver  operating  with  a  linear  channel. 
In  some  practical  applications,  this  may  be  desirable  from  consideration 
of  dynamic  range  requirements.  The  general  expression  for  the  bit  error 
rate  for  both  spread-spectrum  and  PSK  cases  may  be  extended  to  include 
this  situation  by  setting  the  down-link  noise  equal  to  zero. 


The  expression  for  the  bit  error  probability  may  also  be  used  to 
treat  the  case  of  code-division  multiple  access  (CDMA),  involving  a  large 
number  of  spread-spectrum  signals  at  the  limiter  input.  The  noise  source 
at  the  limiter  may  be  regarded  as  the  sum  of  all  the  undesired  signals 
entering  the  limiter  in  addition  to  the  desired  signal.  The  accuracy  of 
the  results  would,  of  course,  be  dependent  upon  how  closely  the  amplitude 
distribution  of  the  sum  of  the  undesirod  signals  at  the  limiter  input 
represented  a  stationary,  zero-mean  Gaussian  distribution.  The  justifica¬ 
tion  is  usually  provided  by  invoking  the  Central  Limit  Theorem  of  proba¬ 
bility  theory  if  the  number  of  signals  at  the  limiter  input  is  large  and 
the  constant  EF  reference  phase  of  each  3ignal  is ^ ndependent  of  all 
others. 
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II  MATHEMATICAL  ANALYSIS 


A.  Error  Probability 

At  the  receiver  output,  the  decision  whether  a  mark  or  space  has 
been  transmitted  is  made  on  the  basis  of  whether  the  filter  output  (see 
Figure  1)  at  the  sampling  instant  is  positive  or  negative.  If  it  is 
assumed  that  the  transmission  of  a  mark  or  a  space  signal  is  equally 
probable  and  that  the  noise  at  the  low-pass  filter  output  has  a  zero-mean 
value,  then  the  probability  of  error,  P  ,  in  a  decision  is  equal  to  the 
error  probability  in  either  a  mark  or  a  space  signal.  Thus, 

P  =  P  =  P 

e  em  es 

The  probability  of  error  in  detecting  a  mark  is  equal  to  the  probability 

that  the  filter  output  will  be  negative  at  the  sampling  instant.  Thus, 

P  is  given  bv 
e 

o 

P  -  P  =  /  p(z)  dz  ,  (1) 

e  em  I  ’ 


whore  p(z)  is  the  probability  density  function  of  the  filter  output.  To 

calculate  P  ,  one  needs  to  determine  p(z),  which  may  be  obtained  by  first 
© 

calculating  the  characteristic  function  of  the  receiver  output,  z  and  then 

taking  the  Fourier  transform.  Alternatively,  P  can  be  obtained  directly 

e 

from  the  characteristic  function  of  Z.  In  the  Appendix  it  is  shown  that, 
if  Z  is  any  random  variable  of  probability  density  function  p(z)  and 
characteristic  function  CMv),  its  cumulative  distribution  function,  P(a), 
is  given  by 


a  00 

(a)  =  j  p(z)  dz  =  \  "  Im  [< 


C  (v)’e 
z 


-iva  dv 


where  I  denotes  that  only  the  imaginary  part  has  to  be  taken.  Corapari- 
m 

son  of  Eqs.  (1)  and  (2)  shows  that 


P  =  P(C) 
e 


00 

=i.i f  x  r 

2  rr  i  m  L 


c2(v, 


Thus,  the  plan  of  approach  for  calculating  P  will  be  to  determine  C  (v) 

e  z 

and  then  to  perform  the  integration.  It  should  be  noted  that  Eq.  (3)  is 

a  general  result  which  can  be  used  to  determine  P  for  a  linear  as  well 

e 

as  a  nonlinear  channel  of  arbitrary  transfer  characteristic. 


B.  Determination  of  C  (v) 

_ z 

To  derive  the  mathematical  expression  for  C  (v),  it  will  be  neces- 

z 

sary  to  obtain  the  expression  for  the  receiver  output  Z.  To  achieve 
this,  we  write  an  expression  for  the  limiter  input  and  then  systematically 
proceed  to  develop  the  expression  for  the  receiver  output. 

The  input  tc  the  limiter  is  assumed  to  consist  of  two  components, 
the  signal  component  and  the  up-link  Gaussian  noiss. 


1.  The  Signal  Component 


The  signal  component  is  represented  as 


s(t)  =  A  cos  [u>  t  +  $(t)  +  0(t ) ] 
o 


During  any  bit  interval,  the  information  modulation,  0(t),  is  either  0 
or  depending  upon  whether  a  mark  or  a  space  is  being  transmitted. 
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Here  $(t)  is  the  pseudo-random  code,  which  enables  the  receiver  to  re¬ 
cover  the  desired  signal  in  the  presence  of  interference  and  noise. 

Since  the  keying  rate  of  the  code  is  generally  several  orders  of  magni¬ 
tude  greater  than  the  information  rate,  the  required  transmission  band¬ 
width  is  usually  very  much  larger  than  is  required  by  information  modula¬ 
tion  alone.  The  analysis  is  independent  of  the  specific  form  of  §(t); 
typically,  the  phase  coding  used  in  digital  systems  is  either  binary  or 
quaternary. 


2.  Up- Link  Gaussian  Noise 


The  expression  for  up-link  narrowband  Gaussian  noise  is 


n(t)  =  x(t)  cos  u)  t  -  y(t)  sin  u>  t  , 
o  o 


where  x(t)  and  y(t)  represent,  respectively,  the  in-phase  and  quadrature 
components  of  the  noise.  For  analysis  it  is  convenient  to  represent  n(t) 
relative  to  the  code  of  the  desired  signal: 


n(t)  =  x  (t)  cos  to)  t  +  $(t)]  -  y  (t)  sin  [tu  t  4-  $(t>]  .  (6> 

1  o  1  o 


No  loss  of  generality  occurs  in  this  procedure,  provided  that  x^t)  and 
y^(t)  are  treated  appropriately.  They  are  independent,  identically  dis¬ 
tributed  Gaussian  processes,  related  to  x(t)  and  y(t)  by  the  following 


expressions: 


x^t)  =  x(t)  cos  $(t)  +  y(t)  sin  $(t) 


y  (t)  =  y (t)  cos  4(t)  -  x(t)  sin  $(t)  .  (8) 


It  is  important  to  note  that  x^t)  and  y^(t)  are,  in  general,  neither 
stationary  nor  independent  of  the  desired  signal,  since 

R^t  -  t')  =  R(t  -  t'>  cos  [*(t>  -  $(t')]  ,  < 


7 


f  (t)  =  A  cos  [in  t  +  $(t)  +  9(t)]  +  x  (t)  cos  [u>  t  +  $(t)] 
i  o  1  o 


-  y  (t)  sin  tu)  t  +  4(t)  ] 
1  o 


=  R(t)  cos  [oj  t  +  $(t)  +  <p(t)] 
o 


(12) 


where  the  envelope,  R(t),  and  the  phase,  cp(t),  are  given  by 


R(t)  =  [A  cos  9(t)  +  x,(t)]2  +  y2(t) 


1 


(13) 


cp( t )  =  arc  tan 


yx(t) 


A  cos  B(t)  +  x^(t) 


(14) 


The  bandpass  limiter  is  assumed  to  be  ideal  in  the  sense  that 


its  output,  f^(t),  is  given  by: 


f  (t)  =  cos  [uj  t  +  $(t)  +  cp(t)] 
o  o 


(15) 


i.e.,  the  envelope  variation  has  been  completely  removed  without  dis¬ 
torting  the  phase  modulation. 


The  signal  is  then  transmitted  to  the  ground  receiver,  and 
noise  is  added  on  the  down  link  and  in  the  receiver  front  end.  Assuming 
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that  this  noise  is  also  stationary,  zero-mean,  and  Gaussian,  the  receiver 
input  may  be  expressed  as 


g(t)  =  a  cos  [u>  t  +  $(t)  +  cp(t)]  +  u<t)  cos  [uj  t  +  $(t)]  -  v(t)  sin  [w  t  +  $(t)' 
o  o  o 

(16) 


where  a  is  a  constant  determined  by  the  amplification  in  the  satellite 
repeater  and  the  losses  on  the  down  link.  Note  that  the  down  link  noise 
is  also  expressed  relative  to  the  code  of  the  desired  signal,  analogous 
to  Eq.  (6). 

The  correlation  operation  at  the  receiver  involves  multiplica¬ 
tion  of  the  receiver  input,  g(t),  by  the  synchronized  receiver  reference 
2  cos  t  +  $(t)],  followed  by  appropriate  low-pass  filtering,  which  is 
equivalent  tc  an  integration  operation  over  tne  bit  interval  T  of  the 
data  signal.  Neglecting  the  filterable  double-frequency  term,  the  re¬ 
ceiver  output  is  given  by: 

T 

Z  =  ~ J"  [a  cos  cp(t)  +  u(t)]  dt  ,  (17) 

o 

where  of  is  a  correlator-gain  normalizing  constant. 

To  avoid  the  mathematical  difficulty  involved  in  solving  the 
above  stochastic  integral,  we  follow  the  approach  used  in  Refs.  1,  2, 
and  3  and  approximate  the  integral  by  a  sura  taken  at  intervals  AT  equal 
to  1/W.  where  W  is  the  bandwidth  of  the  repeater  equal  to  the  reciprocal 
of  the  correlation  time  of  the  p-n  code.  Furthermore,  we  select 
Of  =  /n7a,  which  normalizes  the  receiver  output  noise  power  due  to  up-link 
noise  to  half  in  the  absence  of  signal  and  down-link  noise.  Thus,  Eq.  (17) 


i 


4 


becomes 


vmuxcrmiwt. 


kZ 


COS  ‘J  +  — 

k  a 


*£-\ 


n  =  TW  , 


=  Z  +  Z 

1  2  ' 


where  n  represents  the  TW  product  of  the  systen.  Clearly,  Z  can  be  re¬ 
garded  as  composed  of  the  sum  of  two  statist vcally  independent  variables, 
Z  +  Z, ,  where  Z  is  the  correlator  outrut  due  to  the  limiter  output 
signals,  and  is  the  resulting  output  caused  by  the  receiver  down-link 
noise. 

Since  the  uncorre la noise  processes  at  both  the  limiter  in¬ 
put  and  the  receiver  input  o're  assumed  to  be  Gaussian  and  to  have  a 
symmetrical  power  density  ipectrum  with  respect  to  the  carrier  frequency, 
samples  of  limiter  input  (signal  plus  noise)  and  down-link  receiver  noise 
spaced  1/W  se'  Is  ^;jjtrt  will  be  statistically  independent.  Furthermore, 
the  effect  .mAing  is  merely  to  remove  the  amplitude  modulation 
without  d?_...  ing  the  phase  modulation;  therefore,  statistical  inde¬ 
pendence  will  also  be  maintained  between  the  samples  at  the  limiter 

output .y  Thus,  it  can  be  assumed  that  in  Eq.  (18)  and  u  represent 
f  k  k 

identically  distributed  statistically  independent  random  variables 
spaced  at  intervals  1/W  seconds  apart. 

The  characteristic  function  of  the  sum  of  two  independent 

quantities  is  equal  to  the  product  of  the  individual  characteristic 

function.  Thus,  C  (v)  is  given  by 
z 


C  (v)  =  C  (v)*C  (v) 

z  z  z 

1  2 
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where  C  (v)  and  C  (v)  are,  respectively,  the  characteristic  functions 
Z1  Z2 

of  the  received  signal  and  the  Gaussian  down-link  noise  in  Eq.  (19).  If 
2 

<7  is  the  total  down-link  noise  power  at  the  receiver  input. 


C  (v)  =  e 
Z2 


2  2,2 

-v  a  /2a 


(21) 


In  Eq.  (18)  represents  a  sum  of  n  identically  distributed 

sinusoids.  Therefore,  C  (v)  can  be  written  as 

Z1 


C  (v)  =  [C(v)l 
Z1 


(22) 


where 


2TT  .  V 
C(v)  =  j  e 
o 


‘Tf  ec“  * 


p(cp)  dtp 


(23) 


The  probability  density  function  p(«p)  of  the  phase,  when  a  mark  is  trans¬ 
mitted,  i.e.,  8(t)  =0,  in  Eq.  (14),  is  expressed  by  the  relationship4 


2  2 

-A  /2 a  ®  ,  2  2 

1  r  ~(r  -  2Ar  cos  cp)/2a. 

p(cp)  =  2 - - —  /  re  i 

2TTCT  J 

1  O 


dr 


(24) 


Here  cr^  ropresents  the  total  noise  power  at  the  limiter  input.  Expanding 
the  exponential  function  in  Eq.  (23)  in  a  power  series  yields 


00 

m  / ,  v  P 

— E*»  ■ 


(25) 


P— o 


where 


2tt 

W 


cos*cp  p(cp)  d<p  , 


(26) 


th 

is  the  p  moment  of  ccs  <p. 


11 


12 


These  six  semi-invariants  will  suffice  for  the  expansion  of  C  (v)  through 

-2  Z1 
terras  of  order  n  . 


The  expression  for  C  (v)  in  Eq.  (29)  can  also  be  written  as 

Z1 


C  (v)  =  e 
z 

1 


,  2  cu  oo  _r 

(i  /n  X  v-X  v  /2)  ___  m  ^  X  ,  ,r 
1  2  V  v  n_  j  \  r  /  iv  \ 

ZmJ  mi  rl  \/n / 

m=o  |_r=3 


The  exponential  factor  is  recognized  as  the  characteristic  function  of  a 

Gaussian  random  variable  of  mean  /n  X  and  variance  X  .  The  expression 

1  ,2 

for  C  (v)  ; involves ' the  product  of  the  Gaussian  characteristic  function 
Z1 

and  the  powers  of  v.  It  can  easily  be  shown  by  expanding  Eq.  (31)  that, 

'  i 

following  the  leading  Gaussian  term,  the  successive  higher-order  terms 

1/2  3/2  2 

in  the  series  are  inversely  proportional  to  n  ,  n,  n  ,  n  ,  and  so 


forth.  Thus, 


C  (v)  =  e 
Z1 


(i  /n  X  v-X  v  /2) 
1  ^ 


i  X3  3  1  /  X4  4  X3  6\ 

1-Afrv  +;\iiv  -™v) 


/x  xx 

_i _ I  _5_  y5  _  3  4  y7 

3/2  \ 120  V  "  144  V 


3  9 

+  -  v 

1296 


1  /  X6  6  X4 

2  \  720  V  +  1152 

n  ' 


2  2  4 

4  8  X3X5  8  X4X3  10  X3  12\ 

52  720  1728  31104  I 


C.  Determination  of  the  Moments 


The  six  moments  required  for  the  determination  of  the  semi-invariants 
may  be  obtained  by  first  expanding  cosPcp  in  Eq.  (26)  m  a  Fourier  cosine 
series. 
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p  ev  n 


p 

cos  cp  = 


i— o 


i(P-l) 


P'. 


i=o 


(p  -  i)’.-i’.‘2* 


— r  cos(p  -  2i)  cp  ,  p 
3-1 


odd 


(33) 


and  then  evaluating  the  integral  in  Eq.  (26)  by  using  the  result: 

2n  rl£  _  ji  +  i 


I 


p-2i 


r„  ~(f  - 1  * 1 ) 

cos  (p  -  2i)  cp  p(cp)  dcp  =  _  2i  +  l) 

'  »Fi  (■§  -  ^  p  - 21  *  l>  -p? )  . 


(34) 


where 


2  2,2 

P  =  A  /2a 


(35) 


is  the  limiter  input  SNR;  T(  )  is  the  gamma  function;  and  ^(a,  b>  _x> 
is  the  confluent  hypergeometric  series  defined  as 


2 

a  x a(a  +  1)  >  x_ 

^(a,  b,  -x)  =  1  -  b  j,  +  b(b  +  i)  21 


(36) 


For  even  values  of  p,  the  hypergeometric  function  in  Eq.  (34)  can  be  ex¬ 
pressed  in  closed  form  in  terms  of  x  and  e  X;  while  for  odd  values  of  p 
it  is  expressible  in  terms  of  modified  Bessel  functions,  Iq  and  1^  x  and 
e"^2,  with  the  relationships  given  in  Table  1,  the  following  expressions 
can  be  derived  for  the  six  nionents. 
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The  dependence  of  the  six  moments  and  that  of  the  semi-invariants 

2 

on  the  limiter  input  SNR  are  shown  in  Figures  2  and  3,  respectively. 

These  are  obtained  by  numerically  evaluating  Sqs.  (3?)  and  (30).  For 
2 

large  all  moments  approach  unity,  while  the  semi-invariants,  except 

X^,  which  is  identical  to  m^,  approach  zero.  These  results  are  consistent 

with  those  to  be  expected.  The  contribution  of  the  limiter  input  noise 

to  the  receiver  output  must  decrease  with  an  increase  in  i..put  SNR.  In 
2 

the  limit  when  p^  approaches  infinity,  only  the  signal  component  remains, 

and  there  is  no  contribution  from  the  up-link  noise;  i.e.,  all  X  =0, 

r 

for  r  =  2,  3,  4  . . .  . 

D.  Calculation  of  P 
_ o 

Substitution  of  Eq.  (21)  and  the  imaginary  part  of  Eq.  (32)  in 
Eq.  (3)  yields  the  following  series  for  the  bit-error  probability: 


FIGURE  3  SEMI-INVARIANTS  AS  A  FUNCTION  OF  LIMITER  INPUT  SIGNAL-TO-NOISE- 
POWER  RATIO 
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where 


2  cf 

a  =  —  +  X 
2  2 
a 


(39) 


The  solution  of  the  first  Integral  is  given  in  terns  of  error 
function,7  and  all  the  other  integrals  may  be  evaluated  with  the  aid  of 
the  following  two  integrals: 
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+ - \  H  <p)  + - rH  (p)  -  ...  , 

1728  or  9  31104  a  11 


where  erf(  )  is  the  error  function.  The  Hermlte  polynomials  H  (p)  are 

o 

n 

given  In  Table  2. 

Following  the  leading  error  function  term,  the  first  correction 

term  Is  Inversely  proportional  to  the  square  root  of  the  TW  product  n, 

and  the  successive  higher-order  terms  are  Inversely  proportional  to  n, 

3/2  2 

n  ,  n,  and  so  forth.  Such  an  expansion  is  well  known  as  an  Edgeworth 
series.  The  first  term  in  the  series  yields  the  probability  of  error 
resulting  from  the  receiver  output  component  that  Is  Oaussian-dlstrlbute'd 
the  succeeding  terms  provide  the  contribution  of  the  non-Qaussian  com¬ 
ponents.  The  properties  of  an  Edgeworth  series  have  been  investigated 
in  detail  by  Cramir,*  who  has  shown  that,  under  fairly  general  conditions 
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Table  2 


HERMITE  POLYNOMIALS 


H  (x) 

eo 


H  (x) 
C1 


H  (x)  x  -  1 
e2 


H  (x)  x~  -  3x 
e3 


4  2 

H  (x)  x  -  6x  +  3 
e4 


,  5  3 

H  (x)  x  -  lOx  +  I5x 

e5 


H  (x)  x6  -  15x*  +  45x2  -  15 
e6 


7  5  3 

H  (x)  x  -  21x  +  105x  -  105x 
®7 


8  6  4  2 

H  (:;)  x  -  28x  +  210x  -  420x  +  105 

e 

8 


9  7  5  3 

H  (x)  x  -  36x  +  378x  -  1260x  +  945x 

e9 


10  -8  6  4  2 

H  (x)  x  45x  +  630x  -  3150x  +  4725x  -  945 

e 

10 


11  9  7  5  3 

H  (x)  x  -  55x  +  990x  -  6930x  +  17325x  -  10395x 

G11 


112  in  8  r  4  o 

H  (x)  x  -  66x  +  1485x  -  138(50x  +  51975x  -  62370x"  +  10395 

G  I 


with 
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the  series  gives  an  asymptotic  expansion  of  P  in  powers  of  n 


•1/2 


a  remainder  term  of  the  same  order  as  the  first  term  neglected.  Since 

-2 

we  have  included  the  terms  up  to  the  order  n  ,  the  truncation  error 

may  be  expected  to  be  of  the  order  of  the  next  term,  which  would  be 

-5/2  2  2 

proportional  to  n  .  For  the  range  of  values  of  p  ,  p  ,  and  n  that 

1  2 

would  be  of  interest  in  practical  applications,  it  is  expected  that  the 

first  five  terms  of  the  Edgeworth  series  expansion  will  given  an  accurate 

approximation  of  the  error  rate.  For  most  purposes,  the  expansion  of 

P  through  terms  of  1/n  will  suffice;  however,  for  very  low  error  rates 
6  -5 

(<10  )  the  higher-order  terms  will  be  needed. 


The  contribution  of  higher-order  terms  in  Eq.  (45)  is  primarily 

2 

dependent  on  the  values  of  p^  and  n.  Generally,  when  n  is  small  (between 

20  to  50)  the  up-link  SNR  must  be  high  (around  0  dB)  to  achieve  low 
-5 

error  rates  (<10  ).  Under  these  conditions,  it  is  essential  to  con¬ 

sider  the  higher-order  terms,  since  their  contribution  to  the  error 
rate  will  be  significant.  As  the  processing  gain  n  is  increased  (by 

raising  the  code  chip  rate),  the  required  repeater  bandwidth  W  also 

2 

increases,  and  consequently  the  up-link  SNR  decreases  as  a  result 
of  the  enhancement  in  the  noise  powe~  brought  about  by  the  increased 
repeater  bandwidth.  In  the  limit,  as  n  becomes  very  large,  all  the 
higner-order  terms  in  Eq.  (*.5)  approach  zero,  and  Pg  is  given  by  the 
leading  error- function  term  resulting  from  a  Gaussian  distribution 
of  noise  at  the  receiver  output: 


P  =  7  [1  -  erf(p//2)] 
e  2 


(46) 


2 

For  large  n,  the  up- link  SNR  p  will  also  approach  zero;  Eq.  (43)  thus 
reduces  to 
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(47) 


2 

P 


2  2 

VP2 


1  +  0. 


-» 0 


These  asymptotic  results  for  large  n  are  in  complete  agreement  "’ith  the 
results  obtained  by  Aein.3 

For  a  linear  system  the  error  probability  in  detecting  a  spread- 
spectrum  signal  in  the  presence  of  Gaussian  noise  is  given  exactly  by 
Eq.  (46).  The  corresponding  expression  for  the  receiver  output  SNR  is 


2 

o 


2n 


VP2 


,  2  2 
1  +  P1  +  p2 


(48) 


Thus,  when  n  is  large,  the  distribution  of  the  up-link  noise  at 
the  receiver  output  for  a  hard-limited  system  will  be  approximately 
Gaussian,  provided  also  that  the  up-link  SNR  is  small.  Comparison  of 
Eqs.  (47)  and  (48)  shows  that  under  these  conditions  hard- limiting 
degrades  the  system  performance  by  a  constant  factor  of  tt/4  or  -  1.05 
dB. 

When  the  up-link  SNR  is  high,  the  error  rate  will  be  determined 

primarily  by  the  Gaussian  down-link  noise.  The  higher-order  terms  will 

be  small  compared  to  the  leading  error-function  term,  since  all  the 

2 

semi-invariants  except  X  approach  zero  at  large  values  of  p  ,  as  can 

2 

be  seen  from  Figure  3.  In  the  limit,  as  p^  becomes  very  large,  Pg  is 
given  by  Eq.  (46),  and  Eq.  (43)  reduces  to: 


2  2  2 
P  »2n  p2  ,  pi  . 

The  same  asymptotic  result  is  obtained  for  a  linear  system  [Eq. 
as  wouid  be  expected. 


(49) 

(48)), 
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The  expression  for  the  error  rate  can  also  be  r.3ed  for  the  code¬ 
division  multiple  access  (CDMA)  case,  involving  a  large  number  of  constant- 
envelope  phase-coded  spread-spectrum  carriers  at  the  limiter  input.  The 
noise  source  at  the  limiter  is  then  the  sum  of  all  the  undesired  signals 

entering  the  limiter  in  addition  to  the  desired  signal.  The  up-link 
2 

SNR  is  thus  the  ratio  of  the  desired  signal  power  to  the  sum  of  the 
powers  of  all  the  undesired  signals  and  the  satellite  repeater  noise. 

The  accuracy  of  the  result  is,  of  course,  dependent  on  how  closely  a 
finite  sum  of  p-n  carriers  can  be  modeled  as  a  stationary,  zero-mean 
Gaussian  process. 

For  some  practical  applications,  it  may  be  desirable,  from  con¬ 
sideration  of  dynamic  rarge  requirements,  to  incorporate  a  limiting 

front  end  in  a  receiver  operating  with  a  linear  channel.  The  probability 

2 

of  error  in  this  case  is  obtained  by  considering  a  equal  to  zero,  i.e., 

2 

no  down-link  noise  (o  =  X.  )  in  Eq.  (45)  . 

at 


Ill  DETECTION  OF  A  BIPHASE  PSK  SIGNAL 


The  error  probability  in  the  detection  of  a  hard-limited  binary 
phase- shift-keyed  signal  can  be  determined  by  employing  essentially 
the  same  analytical  approach  as  is  used  for  the  spread-spectrum  signal 

i 

In  this  case,  the  signal  at  the  limiter  input  is  a  constant-envelope, 
sinusoidal  carrier 

i  i 

s(t)  =  A  cos  [u)Qt  +  0(t)]  ,  .  (50) 

where  the  binary  phase  coding  9(t)  is  either  0  or  rr,  depending  upon 
whether  a  mark  or  a  space  is  being  transmitted.  The  above  expression 

i 

corresponds  to  $(t)  =  0  in  Eq.  (4).  If  T  is  the  duration  of  each  data 
bit,  the  information  bandwidth  is  1/2T;  the  required  RF  transmission 

I 

bandwidth  W  is  equal  to  1/T;  and  the  time- bandwidth  product  is  one. 

'  i 

The  calculation  of  the  receiver  low-pass  filter  butput  follows; 

I 

it  is  analogous  to  that  for  the  spread- spectrum  case.  The  filter  out¬ 
put  is  given  by 

! 

u 

Z  =  cos  cp  +  ~  ,  ,  (51) 

a  ,i 

which  corresponds  to  n  equal  to  one  in  Eq.  (13). 

To  determine  the  error  probability  by  using  Eq.  (3),  we  need  the 
characteristic  function  of  the  filter  output.  This  is  given  by 

2  2  2 

_  ,  .  _  r  ivcoscpl  -v  a  /2a  (52) 

c/v>  - E  le  j  • 8 
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where  E[  ]  represents  the  expected  value  or  expectation.  Substituting 
the  imaginary  part  of  Eq.  (52)  in  Eq.  (3)  yields: 


i  b 

■HI 


E[sin(vcos  cp> )  e 


2  2  2 
-v  0  /2a  dv 


(53) 


•l  -  \  eH/S  008  »)]  .  > 


(54) 


!  since  the  integral  in  Eq.  (53)  can  also  be  expressed  in  terms  of  the 
error  function: 7 


00 

•"(ts 008  V)  m  n  /  : 


’  2  2  2 
sin(vcos  ki)  -v  0  /2a  dv 
- x_  e 

V 


(55) 


By  using  the  .expression 


sin(vcoscp)  =  2  \ ^  (-  l)n  •  J  (v)  cos(2n  +  1)  cp  ,  (56) 

j  mtmo^  2m*l , 


n=o 


Eq.  (,55)  becomes 


erf(/t5  C°S  9 


»n/ 


n=o 


J  (v)  22  2 

2n  +1  -vo  /2a  dv-cos(2n  +  1)  cp 

-  e 


(57) 


The  above  integral  is  attributed  to  Weber  and  Sonine;  its  solution  is 
given  in  terms  of  a  confluent  hypergeometric  series.9  The  result  is: 


i 
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2n  +  2, 


cos  (2n  +  1)  cp  •  p(cp)  dcp 


(59) 


The  expected  value  of  the  cosine  function  in  Eq.  (59)  can  be  evaluated 
with  the  aid  of  Eq.  (34): 

2tt 

cos  (2n  +  1)  cp  p(cp)  dcp 
o 

(60) 


Substituting  Eq.  (60)  in  Eq.  (59),  we  obtain  the  following  expression 
for  the  probability  of  error: 
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2  TT 


(-l)n-r(n  +  M-r(n  +  -) 

\  2/  \  2/  2n+l 

- L. - L  p 

2  K1 

[ (2n  +  1)!] 


1F1  (”  +  i  >  2n  +  2>  -Pi)  •  !Fi  (»  +  1  -  2"  +  2>  -Pa)  -  . (61) 


where 


2  2  2 

=  A  /2o^  =  limiter  input  SNR, 


2  2  2  2 

=  a  /2o  =  P^/c  =  down- link  received  satel lite- power- to-noise- 

power  ratio. 

Here  ?  is  the  satellite  output  power  referenced  to  the  ground  terminal 
and  contains  both  signal  and  up-link  noise  as  delivered  by  the  limiter. 
The  series  in  Eq.  (61)  is  corvergent,  which  may  be  shown  as  follows. 

We  have  the  general  relation 


|jF^  (a,  b,  -x)  j  £  i  ,  Re  a  >  0  . 


Hence,  the  series  in  Eq.  (61)  will  converge,  if  the  dominating  series 

obtained  by  setting  both  the  confluent  hypergeometric  functions  equal 

to  one  converges.  The  resulting  series  is  absolutely  convergent  for 
2  2 

all  values  of  p  and  p  .  This  can  easily  be  shown  by  D'Alembert's 
1.  2 

ratio  test  for  absolute  convergence: 


n-*«>  U 


=  Lim 


l\(  3  \ 

n  +  i)\”  +  i)  f  ,2 

~~T~~~~2  lpip2>  -0 


n-*°  2 

(2n  +  3)  • (2n  +  2) 


1  if,  ;,-l'f'  k  1  '  ■ '  -. - —  -  ■■■'■.*  ...  i-b  A'_  i  .  i.*  ■  ■*  ■  .’*:•>■  ,.  re,..->  O-'T-  . 


An  alternative  expression  for  P  which  is  more  suitable  for  numeri- 

e 

cal  computation  is  obtained  by  replacing  the  confluent  hypergeometric 
function  by  the  relationship 


1F1  (n  +  \  >  2n  +  2' 


-  P 


2\  e 


•p2/2 


*  n;2 


_2n 


2n 


[ln  ( p2/2)  +  I^tp2^)]  .  (64) 


The  exnression  for  P  then  becomes: 

e 


i  V’a  ‘(V^)72 


V1  (-*>r 

jLj  (2n 


n=o 


1) 


>„®  *  ■  J) 


(pj\  (>l) 

1A2  )+  Iftfi\2  / 


(65) 


It  is  interesting  to  examine  the  behavior  of  P  for  large  SNRs. 

e 

2 

When  is  large,  the  confluent  hypergeometric,  function  in  Eq.  (61) 
may  be  replaced  by  its  asymptotic  expansion  for  large  negative  argument. 
The  result  is: 


,  -  (-D 

1  _ 


P  = - -  > 

e  2  tt  /  -4 


°T(n  +  i) 


2n+l 


(2n  +  1).’  y2 


iTi  (” + i 2n  +  2>  -  4) 


n=o 


=  -  [1  -  erf  p2]  ,  pA  ->  »  , 


(66) 
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since  the  series  represents  1/2  erf  p  ,  as  may  be  seen  from  Eq.  (58). 

Equation  (66)  is  the  expected  result  for  ideal  limiting  in  the  absence 

2 

of  up-link  noise.  Similarly,  when  p  is  large,  P  is  given  by: 

2  e 


P 

e 


=  2  (1  "  erf  P]>» 


0  ->>  co 


(67) 


This  is  identical  with  the  expression  for  the  error  probability  in 
detecting  a  PSK  signal  in  the  presence  of  Gaussian  noise  by  using  a 
correlation  detector.  Thus,  incorporation  of  a  limiting  front  end  in 
a  receiver  operating  with  a  linear  channel  will  not  degrade  the  signal 
detectability  compared  to  a  linear  correlation  receiver. 
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IV  NUMERICAL  RESULTS 


A.  Output  SNR 

2 

Figures  4  and  5  show  the  receiver  output  SNR  p  [Eq.  (43)]  as  a 

2  2 

function  of  up-link  and  down-link  SNR  p  and  pr,  respectively.  The 
dashed  lines  are  similar  curves  for  a  linear  system.  The  output  SNR 
for  a  linear  system  is  given  by  Eq.  (48). 

2 

At  low  values  of  p  ,  Eq.  (43)  reduces  to  Eq.  (47).  Comparison  of 

Eqs.  (47)  and  (48)  shows  that  the  presence  of  a  limiter  degrades  the 

output  SNR  by  a  factor  of  tt/4  or  -1.05  dB  compared  to  a  linear  system. 

2  2 

As  increases,  the  degradation  in  p  decreases,  and  eventually  there 

is  an  improvement  in  the  receiver  output  SNR  over  that  of  a  linear  sys- 

2 

tem  (Figure  4)  .  jut*t  p  ‘  represent  the  up-link  SNR  that  will  yield  the 
2 

same  value  for  p  for  both  a  hard-limited  and  a  linear  system  for  a 
2 

given  value  of  p  . 

2  2 

Figui  4  shows  that  for  p^  >  p^there  is  an  improvement  in  the  SNR 

2 

for  a  hard-limited  system.  However,  as  p^  becomes  larger,  the  improve- 

2 

ment  gets  smaller.  This  is  because  p  is  now  determined  by  the  down¬ 
link  noise.  The  same  output  SNR  is  obtained  for  both  a  hard-limited 
and  a  linear  system  when  there  is  no  up-link  noise.  Figure  4  also  shows 

that  p  increases  with  decreasing  p".  The  minimum  value  (approximately 
10  2 

-4  dB)  is  obtained  when  there  is  no  down- link  noise. 

In  the  absence  of  down-link  noise  and  a  large  up-link  SNR,  Eq.  (43) 
simplifies  to: 

2  4  2  2 

P  as Snp,  ,  P2  =  ^  >  Px  00 
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LIMITER  INPUT  SIGNAL-TO-NOISE  RATIO  P  2  —  dB 
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FIGURE  4  RECEIVER  OUTPUT  SIGNAL-T0-N0I3E-P0WER  RATIO  AS  A  FUNCTION  OF 
LIMITER  INPUT  SIGNAL-TO-NOISE-POWER  RATIO 


IN  dB  -10  log 
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FIGURE  5  RECEIVER  OUTPUT  SIGNAL-TO-NOISE-POWER  RATIO  AS  A  FUNCTION  OF 
DOWN-LINK  RECEIVED  SATELLITE-POWER-TO-NOISE-POWER  RATIO 
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Davenport  has  shown  that  a  bandpass  limiter  degrades  the  input  SNR, 
2  2 

p  ,  by  a  factor  of  tt/4  at  low  values  of  but  provides  a  constant 


improvement  of  3  dB  at  high  values  of  p^.  Equation  (43)  shows  that,  if 


the  bandpass  limiter  is  followed  by  a  correlation  detector,  the  degrada¬ 


tion  in  the  input  SNR  will  remain  tt/4  at  low  values  of  p  ;  however,  at 

2  2 
large  values  of  p^  there  will  be  an  improvement  of  p  ,  which  is  propor¬ 


tional  to  the  square  of  the  input  SNR. 


B.  Probability  of  Error  for  a  Spread-Spectrum  Signal 


The  expression  for  the  probability  of  error  in  the  detection  of  a 


constant-envelope  spread-spectrum  signal  after  transmission  through  a 


hard  limiter  is  given  by  Eq.  (45).  The  semi-invariants  are  functions 


of  up-link  SNR  p  and  are  shown  plotted  in  Figure  3.  Figure  6  shews 
^2  2 

P  as  a  function  of  p  for  constant  values  of  p  .  The  dashed  lines  are 
el  3 


similar  curves  for  a  linear  system.  The  probability  of  error  for  a 


linear  system  is  given  by  Eq.  (46),  where  the  receiver  output  SNR  is 


determined  by  Eq.  (48).  Figure  6  shows  that  P  decreases  monotonically 

6  2 


with  increasing  up-link  SNR;  however,  at  large  values  of  p^  the  perfor¬ 


mance  of  the  system  exhibits  an  irreducible  error  probability  represented 


by  the  bottoming  of  the  error  rate,  which  depends  on  the  SNR  on  the 
down  link.  This  means  that,  as  more  power  is  placed  in  the  information¬ 


bearing  signal,  the  ultimate  performance  of  the  system  is  governed  by 

2 

the  down-link  SNR  For  large  values  of  p^  all  tne  semi-invariants 


except  approach  zero,  and  Eq.  (45)  reduces  to  Eq.  (46),  which  is 


then  identical  with  the  expression  of  P^  obtained  for  a  linear  system. 


Figure  6  shows  that  the  error  rate  for  a  hard-limited  system  is  higher 


than  that  for  a  linear  system  at  low  values  of  p  but  that  it  gets 

2 

smaller  as  the  up-link  SNR  is  increased.  Of  course,  if  p^  became  very 


large,  the  error  rate  for  a  hard-limited  system  and  that  for  a  linear 
system  will  approach  the  same  limiting  value,  which  is  determined  by 


1 


FIGURE  6  ERROR  RATE  AS  A  FUNCTION  OF  UP-LINK  SIGNAL-TO-NOISE-POWER 

RATIO  FOR  CONSTANT  VALUES  OF  RECEIVED  SATELLITE-POWER-TO-NOISE- 
POWER  RATIO 


the  down-link  SNR.  The  improvement  in  the  error  rate  performance  of  a 

hard-limited  system  results  from  the  fact  that  the  receiver  output  SNR 

2  2 
p  is  higher  than  that  for  a  linear  channel,  when  the  up-link  SNR  p 

2  2  * 
has  become  sufficiently  large  (p^  >  p^.  Examination  of  Eq.  (45)  shows 

that  in  this  region  the  major  contribution  to  the  error  rate  is  provided 

by  the  leading  Gaussian  term.  The  successive  higher-order  terms  tend 

to  increase  the  error  rate,  but  their  contribution  is  so  small  that  P 

still  remains  lower  than  for  a  linear  system. 

In  Figure  7  the  error  rate  is  shown  as  a  function  of  the  down-link 

2 

received  satellite-power- to- noise- power  ratio  p  for  different  values 

2  2 
of  p  .  The  dashed  lines  represent  the  error  rate  that  is  obtained  if 

only  the  leading  Gaussian  term  in  Eq.  (45)  is  considered,  nd  the  output 
2 

SNR  p  is  determined  from  Eq.  (43).  The  contribution  of  ,e  higher- 
order  terms  is  the  difference  between  the  solid  and  the  dashed  lines. 

It  can  be  seen  that  this  contribution  would  be  significant  at  low  error 

—5  2  2 
rates  (<10  )  and  moderate  values  of  p  (>0  dB)  .  As  p  increases, 

2  2 

the  error  rate  performance  of  the  system  is  determined  primarily  by 
the  contribution  of  the  up-link  noise  at  the  receiver  output,  and  the 
higher-order  terms  in  Eq.  (45)  must  be  considered.  The  curves  of  Figure 
7  also  show  a  bottoming  of  the  error  .rate,  representing  the  irreducible 
er:’or  probability  brought  about  by  the  presence  of  the  up-link  noise. 

Figure  8  shows  the  error  rate  as  a  function  of  the  processing  gain 

2  2  2 

for  constant  values  of  p  and  two  values  of  p  The  curves  for  p 

1  2  2 

i.e.,  no  down-link  noise,  represent  the  case  of  a  limiting  front  end 
incorporated  in  a  receiver  operating  with  a  linear  channel.  When  the 

_  3 

error  rate  is  lov;  (<10  ),  the  higher-order  terms  in  Eq.  (45)  must  be 

considered,  particularly  when  the  processing  gain  is  small,  and  the  up¬ 
link  SNR  is  greater  than  -10  dB.  This  can  be  seen  clearly  from  Figure 

9,  where  P  is  plotted  as  a  function  of  the  up-link  SNR  for  constant 
e 
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FIGURE  7  ERROR  RATE  AS  A  FUNCTION  OF  RECEIVED  SATELLITE-P0WER-70-N0ISE 
POWER  RATIO  FOR  CONSTANT  VALUES  OF  LIMITER  INPUT  SIGNAL-TO- 
NOISE-POWER  RATIO 
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LIMITER  INPUT  SIGNAL-TO-NOISE  RATIO  P  j  —  dB 

SA-1328-B 

FIGURE  9  ERROR  RATE  AS  A  FUNCTION  OF  LIMITER  INPJT  SIG^AL-TO-NOISE-POvVER 
RATIO  FOR  CONSTANT  VALUES  OF  PROCESSING  GAIN 


l 


S 


values  of  n  in  the  absence  of  down-link  noise.  The  dashed  lines  represent 

the  error  probability  if  only  the  leading  errcr-function  term  in  Eq.  (45) 

is  considered.  If  the  up-link  SNR  is  kept  constant  (by  either  increasing 

the  transmitted  signal  power  or  reducing  the  data  rate),  while  the 

processing  gain  is  steadily  increased,  the  error  rate  decreases  rapidly 

but  the  contribution  of  the  higher-order  terms  becomes  increasingl*' 

significant.  This  may  be  explained  as  follows:  As  the  process,  z  in 

2 

is  increased,  the  output  SNR  p  increases  linearly  with  n,  as  s 

;  ■  i  2  1 

Eq.  (43).  For  p  >6  dB,  the  leading  term  in  Eq.  (45)  may  be  ap. 
mated  by  the  relationship 


i 

whicl  is  obtained  by  replacing  the  error  function  by  the  first  term  of 
its  asymptotic  expansion  for  large  arguments.  Thus,  Eq.  (45>  becomes 


-p2/2 

V  2tt  p 


+  higher- or dev  terms 


An  n  is  increased,  the  first  higher-border  term  increases  linearly  with 
n,  since  p  =  All  the  semi-invariants,  however,  remain  un¬ 

changed,  since  the  up-link  SNR  is  kept  constant.  Similarly,  it  can  be 

shown  that  the  other  higher-order  terms  also  increase.  The  same  will 

2  2 

also  be  true  in  the  presence  of  down-link  noise  if  both  p  and  p  are 

JL  M 

held  constant  while  n  is  increased.  Thus,  under  these  conditions  the 
receiver  output  noise  will  not  approach  a  Gaussian  distribution  even 
at  large  prqcessing  gains,  and  it  will  be  essential  to  consider  the  con¬ 
tributions  of  the  higher-order  terms  in  the  calculation  of  the  error  rate. 
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Investigation  of  Eq.  (45)  has  shown  that  the  series  expansion  up 
-2 

to  the  order  n  gives  an  excellent  approximation  of  P  ,  even  at  very 

*■  5  ~6  ^ 

low  error  rates  (<10  ) .  Only  for  P  <  10  and  n  4  25  will  more  terms 

e 

be  required  in  Eq.  (45).  In  practical  systems  the  processing  gain  is 

usually  larger  (n  >50),  and  Eq.  (45)  may  be  used  to  determine  P  for 

2  2  6 

virtually  all  values  of  p  and  p  . 

X  2 

An  idea  of  tne  effect  of  the  higher-order  terms  in  the  calculation 
of  the  error  rate  may  be  obtained  from  a  few  selected  cases  considered 
in  Table  3.  If  only  the  leading  error-function  terra  in  Eq.  (45)  is 
considered,  the  resulting  error  rate  is  as  shown  under  the  column  en¬ 
titled  One  Term.  The  other  columns  in  Table  3  show  the  modification  of 

P  as  higher-order  terms  are  included.  Thus,  the  last  column  represents 
e 

the  error  rate  if  all  the  terms  in  Eq.  (45)  are  considered.  It  isi  clear 

that  the  higher-order  terms  provide  improved  accuracy  and  convergence 

of  P  ,  particularly  when  n  is  small.  Even  when  n  is  large,  inclusion 
6  -1 

of  the  terms  through  the  order  n  is  desirable  when  dealing  with  low 
error  rates. 

If  the  up-link  transmitted  signal  power  is  held  constant  while 

the  processing  gain  is  steadily  increased  by  raising  the  chip  rate, 

2 

p^  will  decrease  as  a  result  of  the  increase  in  the  noise  power  at  the 

limiter  input,  brought  about  by  the  increase  in  the  transmission  band- 

2 

width  W.  in  the  limit,  as  n  becomes  very  large,  p^  approaches  zero, 

and  Eq.  (43)  reduces  to  Eq.  (47).  In  the  absence  of  down-link  noise, 

2 

p  reduces  to: 

2 

p2  =  2WT  4  *  2^  =  4  *  E/11  ’  (70> 

2 

where  E  =  A  •  T  represents  the  signal  energy.  This  shows  that  an  in¬ 
crease  in  the  processing  gain  will  not  affect  the  output  SNR  and  that 


42 


3377  x  10  2.5062  x  10  2.4861  x  10  2.4847  x  10  2.4848  x  10-  -17  3  1000 


therefore  ail  the  terms  in  Eq.  (45)  containing  p  will  also  remain  un¬ 
affected.  However,  the  effect  of  the  non-Gaussian  terms  will  decrease, 

-1/2  -1 

since  they  are  proportional  to  n  ,  n  ,  and  so  forth  In  the  limit, 
as  n  becomes  very  large,  the  contribution  of  the  up-iink  noise  at  the 
receiver  output  will  tend  to  be  a  Gaussian  distribution,  and  the  error 
rate  will  be  given  by  the  leading  term  in  £q.  (45).  The  presence  of 
down-link  noise  will  not  change  the  Gaussian  distribution  at  the  receiver 
output,  and  the  error  rate  will  still  be  determined  by  the  leading  term 
in  Eq.  (45),  with  the  corresponding  receiver  output  SNR  given  by  Eq. 


(47)  . 


Figures  10  through  13  provide  additional  numerical  results  for  the 


error  rate. 


C.  Probability  of  Error  for  a  PSK  Signa] 

Numerical  evaluation  of  the  error  rate  for  a  PSK  signal  as  a  func- 
2  2 

tion  of  p  ,  with  p  as  the  parameter,  is  shown  in  Figure  14.  Similar 
12  2  2 

curves  for  P^  as  a  function  of  p  ,  with  p^  as  the  parameter,  can  be 

2  o 

obtained  by  merely  interchanging  p  and  p  in  Figure  le,  since  Eqs.  (61) 

2  *2  * 

and  (65)  are  symmetrical  in  p  and  p  .  Both  Eq.  (61)  and  Eq,  (65)  were 

1  M 

programmed,  and  the  results  obtained  were  identical  with  those  expected. 

2  2 

However,  Eq.  (65)  converges  faster  when  both  P  and  p^  are  large  (>6  dB), 
which  /ill  be  the  case  in  practice.  The  dashed  lines  show  the  error 
rate  for  a  linear  system  [Eq.  (46):  n  =  1] .  The  system  exhibits  a 


bottoming  of  the  error  rate,  representing  an  irreducible  error  prob¬ 
ability,  which  depends  on  the  noise  present  on  either  the  up  link  or  the 

2  2 

down  link,  depending  upon  whether  the  abscissa  is  p  or  p",  It  can  be 

2  2  i 

seen  in  Figuie  14  that,  as  p  tends  to  infinity,  both  a  hard-limited 
system  and  a  linear  system  tend  to  the  limit  given  by  Eq.  i.66).  The 
significant  difference,  however,  xs  that  a  hard-limited  system  approaches 
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FIGURE  10  ERROR  RATE  AS  A  FUNCTION  OF  RECEIVED  SATELLIT£-POWER-TO-DOVVN- 
LINK-NOISE-POWER  RATIO  FOR  CONSTANT  VALUES  OF  PROCESSING  GAIN. 
P,2  «  -13  dB 
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RECEIVED  SATEU.ITE-POWER-TO-DOWN-UNK-NOISE-POWER  RATIO  p*  _  dB 
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ERROR  RATE  AS  A  FUNCTION  OF  RECEIVED  SATELLITE-POWER-TO-DOWN- 
LINK-NOISE -POWER  RATIO  FOR  CONSTANT  VALUES  OF  PROCESSING  GAIN, 


FIGURE  13  ERROR  RATE  AS  A  FUNCTION  OF  RECEIVED  SATELLITE-POWER-TO-DOWN- 
IINK-NOISE-POWER  RATIO  FOR  CONSTANT  VALUES  OF  PROCESSING  GAIN, 


-2«  -18  -12  S  0  6  12  18  2« 

LIMITER  INPUT  SIGNAL-TO-NOISE  RATIO  Pf  —  dB 

SA-132B-13 


FIGURE  14  PSK  ERROR  RATE  AS  A  FUNCTION  OF  LIMITER  INPUT  SIGNAL-TO-NOISE 
POWER  RATIO  FOR  CONSTANT  VALUES  OF  RECEIVED  SATELLITE-POWER 
TO-NOISE-POWER  RATIO 
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the  limit  much  faster,  i.e.,  at  lower  values  of  p  than  a  linear  system 

2  2  2 

for  any  constant  value  of  p  .  As  both  p  and  p  approach  infinity, 

Ct  i  M 

i.e.,  perfect  phase  measurement  at  the  receiver,  approaches  zero. 


V  CONCLUSIONS 


The  principal  conclusion  of  this  research  effort  is  that  the  error 

rate  in  the  detection  of  a  constant-envelope  spread-spectrum  signal  after 

transmission  through  a  hard  limiter  can  be  calculated  very  accurately 

by  using  an  Edgeworth  series  expansion.  The  series  provides  an  asymptotic 

-1/2 

expansion  of  the  error  rate  in  powers  of  n  ,  where  n  is  equal  to  the 

TW  product  or  the  processing  gain  of  the  system.  Inclusion  of  the  terms 
-2 

up  to  the  order  n  should  be  fully  adequate  to  calculate  the  error  rate 
for  virtually  all  ranges  of  values  of  up-link  and  down-link  SNRs  and 
for  the  system  processing  gain  normally  encountered  in  practical  applica¬ 
tions.  Even  when  n  is  small  (n  =  25),  Eq.  (45)  containing  the  first 

five  terms  of  the  Edgeworth  series  can  be  used  to  calculate  very  low 

—6 

error  rates  (P^  a  10  ).  As  n  increases,  the  validity  of  Eq.  (45) 

extends  to  even  lower  error  rates.  In  the  limit,  as  n  becomes  very 
large,  Pg  is  determined  by  the  leading  error-function  term  of  Eq.  (45), 
which  is  also  identical  with  the  result  obtained  by  Aein.3 

The  expression  for  the  error  probability  can  also  be  used  for  the 
code-division  multiple  access  (CDMA)  case,  involving  a  large  number  of 
constant-envelope,  phase-coded  spread-spectrum  carriers  at  the  limiter 
input.  The  noise  source  at  the  limiter  is  then  the  sum  of  all  the  un¬ 
desired  signals  entering  the  limiter  in  addition  to  the  desired  signal. 

The  accuracy  of  the  results  is,  of  course,  dependent  on  how  closely 
the  amplitude  distribution  of  the  sum  of  the  undesired  signals  at  the 
limiter  input  represents  a  stationary,  zero-mean  Gaussian  distribution. 

The  justification  is  usually  provided  bv  invoking  the  Central  Limit 

Theorem,  if  the  number  of  signals  at  the  limiter  input  is  large,  and 

the  constant  RF  reference  phase  of  each  signal  is  independent  of  all  othi  rs. 
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The  error  rate  in  the  detection  of  a  PSK  signal  may  be  calculated 

from  either  Eq.  (61)  or  Eq.  (65)  for  arbitrary  values  of  up-link  and 

down- link  SNRs.  However,  Eq.  (65)  appears  to  be  more  suitable  for 

2  2 

numerical  computation,  particularly  when  p  and  p  are  greater  than  6  dB. 

X  Ct 

This  would  normally  be  the  case  in  practice  to  achieve  error  rates  less 

-p 

than  10  ‘  .  It  should  be  noted  that  the  model  for  the  PSK  case  assumes 
that  the  repeater  bandwidth  is  just  wide  enough  to  pass  the  signal  with 
negligible  distortion  and  to  limit  the  input  noise  to  the  bandwidth 
of  the  signal.  This  would  be  the  situation  in  a  channelized  satellite 
repeater,  where  each  channel  was  used  to  transmit  and  limit  a  single 
PSK  signal. 

The  conclusions  reached  in  this  report  are  equally  valid  for  the 
case  of  no  down-link  noise  which  would  represent  incorporation  of  a 
limiting  front  end  in  a  receiver  operating  with  a  linear  channel.  In 
some  practical  applications,  this  may  be  desirable  from  consideration 
of  dynamic  range  requirements. 


Appendix 
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where  C(v)  =  C(-v)  represents  the  complex  conjugate  of  C(v).  Further- 

i 

more,  the  two  functions  in  Eq.  (A-2)  ^re  also  complex  conjugates  of 
each  other.  Thus, 


-ivz  — -  —  ivz 
e  -  C(v)  e 


=  2i  I 


(A-3) 


where  I  denotes  that  the  imaginary  part  must  be  taken. 

m  ( 

Finally,  with  the  aid  of  Eqs.  (A-2)  and  (A-3).,  P(z)  can  be  expressed 

as 

i 
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